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I. INTRODUCTION 



Recently, the CLEO collaboration reported the observation of some nonleptonic B meson 
decays into two light pseudoscalar mesons with the following branching ratios JTIJ 
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And the upper limits of the branching ratios for other decay modes were of the order of 
10~ 5 . These measurements at CLEO and the future experiments at B-factories motivate a 
great deal of theoretical interests in nonleptonic B decays into two light mesons |3|-[T0 . 



The standard theoretical framework to study non-leptonic B decays is based on the ef- 
fective Hamiltonian approach. The short and long distance QCD effects in the non-leptonic 



weak decays are separated by means of the operator product expansion |TT| . The resulting 
effective Hamiltonian consists of products of scale-dependent Wilson coefficients Cj(//) and 
the local four-quark operators at the renormalization scale \x. The short distance contri- 
butions of the coefficients C;(/x) has been evaluated up to the next-to- leading logarithmic ( 



NLL ) order [12,13] . In the NLL precision, the Wilson coefficients depend on the renormal- 
ization scheme as well as renormalization scale. These dependences should be canceled by 
corresponding scheme/scale dependence of the matrix element of the operators. However, 
hadronic matrix elements are usually calculated under the factorization approximation in 
which they are replaced by the scale/scheme independent form factors and decay constants. 
In order to achieve the scale/scheme independence of the matrix elements of the effective 
Hamiltonian, Ali and Greub Q included the corrections of the one-loop penguin-like di- 
agrams and some process independent part of the vertex corrections associated with the 
four-fermi operators to the partonic matrix elements before doing factorization step. They 
also have taken into account the effects of the 0(a s ) tree-level matrix element associated 
with the chromomagnetic dipole operator. Such corrections compensate the scale/scheme 
dependences of the Wilson coefficients and one finally get the effective Hamiltonian involving 
the scale/scheme independent Wilson coefficients and four-quark operator. 

The factorization method is successful to describe the heavy-to-heavy nonleptonic decays 



14j| . The intuitive argument for the factorization is given by Bjorken [fL5[] based on the idea 
of color transparency. The quark pairs with high energy coming from the heavy meson 
decay hadronize after they have traveled some distance from each other. Hence the decay 
process is expected to factorize into the color singlet current pairs since soft gluon effects 
are small after the hadronization occurs. There is more theoretical argument about the 
factorization based on QCD. In the specific kinematic region in which the two light quarks 
are highly collinear and all quarks are almost on-shell, the leading term of the Green function 
in expansion of inverse powers of the heavy quark masses and the large energy transferred 
to the light quark pairs exhibits factorization |I6| . The authors in Ref. Q extended this 
framework to the heavy-to-light transitions. They also introduce a new free parameter £ 
that describes nonfactorization effects such as color octet contributions and concluded that 
the range < £ < 0.5 is consistent with data H. However one should keep in mind that the 
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factorization for the heavy-to-light transition is still only a model to describe the complex 
hadronic matrix elements since there is no theoretical justification as the heavy-to-heavy 



transition case | I6| . In principle there is no reason that only one single value of the parameter 
£ can explain the branching ratios of all kind of different modes. For example, the authors of 
Ref. introduce two different £'s corresponding to the different currents structures which 
give different nonfactorization corrections. However, we will use the factorization method 
including one nonfactorization parameter £ in this work. 

On the other hand, there are some possible anomalies in the inclusive B decays. The 
persistent discrepancy of the measurements of the semileptonic branching ratio and charm 
multiplicity with the theoretical predictions have been known quite for a while [|I7]]. One can 
argue that these problems arise because of the breakdown of local quark-hadron duality that 
was invoked when one estimates the nonleptonic B decay rate. However, several authors 
have noticed that these puzzles can be solved if one assumes that the Wilson coefficient of the 
chromomagnetic operator is enhanced by some new physics contributions jnj. Incidentally, 
this may help to understand the recently measured branching ratio for B(B — > r]'X s ), which 
is surprisingly larger than previously expected. Such an enhanced b —>■ sg should affect 
inevitably the exclusive decay rates of B mesons, and it is one of the main themes of our 
present work to study such effects. 

In this work, we consider the nonleptonic B decays into two light pseudoscalar mesons and 
investigate the possible effects through the enhanced chromomagnetic dipole contribution. 
The branching ratios are calculated by the generalized factorization method proposed in 
Ref. ||. The QCD and EW penguin contributions are included in this work. The two-angle 
mixing formalisms are used in the calculation of the decay rates involving rj, 7/ mesons, and 
the amplitudes b — > s(gg) — > 5(77,77') |J are included using the QCD anomaly, instead of 
the intrinsic charm contents of rj^'\ We also study the CP asymmetry in B meson decays, 
both the direct CP asymmetries of charged B meson decays and the time integrated CP 
asymmetries of neutral B meson decays. In this calculation, we assume that the strong 
phases are given by the penguin-type diagrams with internal light quarks which could have 
on-shell momentum 



II. CALCULATIONAL FRAMEWORK 

The AB = 1 effective Hamiltonian is given by 



H eff - -7= E 

V z q =d,s 



v uh v: q {c x o\ + c 2 o u 2 ) + VcV^c.oi + c 2 o c 2 ) - v tb v; q (e^o, + <3a) 

(2) 

where the operators are 
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7 = l(q a b a ) V -AE q >e q > (q^q p ) v+A 8 = |(g /3 & Q )v-A E,' z<f (C?i) y+A 

o 9 = UvMv-AE q >e q > {q p q p ) v _ A o w = |(# Q )^E 9 'V {qJp) v _ A 

O g = (g s /8n 2 ) m b s a a^(l + 75) 2)fy 

Here the Oi and 2 are the current-current operators. The 3 ~ 6 and 7 ~ Oi are called 
the gluonic and electroweak penguin operators respectively, and O g is the gluonic dipole 
moment operator. The subscripts V ±A represent the projection operators 1 ±75 onto left- 
and right-handed spinor and e q i indicates the electromagnetic charge of the corresponding 
quarks. The subindex a and (3 are the SU(3) color indices and Afg (A = 1 ~ 8) are the 
Gell-Mann matrices. 

If we take mg e = 175 GeV, a s (M z ) = 0.118, a(M z ) = 1/128, we have the following 
numerical values of the Wilson coefficients at the renormalization scale /i = 2.5 GeV in the 
naive dimensional renormalization scheme Q: 

d = 1.117 C 2 = -0.257 C 3 = 0.017 C 4 = -0.044 
C 5 = 0.011 C 6 = -0.056 C 7 = -1 x 10- 5 C 8 = 5 x 10" 4 (4) 
= -0.010 C10 = 0.002 C g = -0.158 

where the Wilson coefficients C\ ~ Cq are taken as the NLL values with respect to QCD. 
The values of the remaining coefficients are given at the leading logarithmic precision. 

In the NLL precision, the matrix elements of local 4-fermion operators Oj's are to be 
treated at the one-loop level. In Ref. ||, the contributions arising from the penguin- type 
diagram of the operators 0\ ~ 0% and the tree-level diagram of the dipole operator O g 
have been calculated and absorbed into the effective Wilson coefficients C^. The process- 
independent contributions from the vertex-type diagrams are also considered. These full 
NLL considerations are sufficient to compensate the scale and scheme dependence arising 
from the factorization approximation and the resulting effective Wilson coefficients are given 
in Ref. |,§] . 

In this paper, we use the factorization approximation in order to calculate the hadronic 
matrix elements of the type < h\h2\Oi\B >, where is a light pseudoscalar meson. In this 
approximation, the hadronic matrix elements are factorized into a product of two matrix 
elements of quark bilinear operators. As an example, let us consider the decay B° — > 7r + 7r~, 
whose matrix element is given by 



M=^{k 6 K>i- V tb V t * d 



2ml(a 6 + a s ) 

O4 + O10 + 



where 



m u + m d )(m b - m u ) _ 
x < 7T~\du-\0 >< n + \ub-\B° > (5) 



Cf f + f C-l({i = odd) , a, = Ctli + i Ct ss {% = even). (6) 
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Here, instead of 1/N C , we introduce a free parameter £ which is supposed to describe the 
nonfactorization contribution. (The amplitudes for other decays are explicitly listed in 
Appendix.) 

When we express the matrix elements of quark bilinear operators in terms of meson decay 
constants and form factors, it is important to remember which convention we use for the 
meson wave function in terms of quark flavor contents. We adopt the following convention, 
for which the meson state is the same as the isospin eigenstate without extra sign. 

\tt + >= \ud>, \n >= \-^=(dd — uu) >, \n~ >= \ud > 

v2 

\K~ >=\-su>, \K°>=\sd>, \K°>=\ds>, \B~ >=\-bu>, \B° >= \bd > 

1 - 1 

\% >= \—=(2ss — uu — dd) >, 1 770 >= |— p(— ss — uu — dd) > (7) 
V6 V3 

With this convention, we define the decay constants and form factors as follows: 

< n~(p)\du_\0 >= %UV^ < K~(p)\su-\0 >= if K P^ 

< Vs(p)\-^(uu^ + dd_ - 2ss_)|0 >= ifep^ 

< r) (p)\-^=(uu- + dd^ + ss_)|0 >= if p^ 

< ri { '\p)\uu-\0 >= if%,)Pn, <V { '\p)\ss-\0>=if° i , ) p l ,, < ^'\p)\cc.\0 >= iffap,* (8) 
and 

< *-(j/)\db.\B-(p) >= l(p + p% - ^^q»]Fr*(<l 2 ) + ^^q.Fr^q 2 ) 

< K-(p>)\sb-\B-(p) >= l(p + pX - ^^ q ,]Fr K (q 2 ) + ^^i^V) 



qz f qz 
2 ™2 ^,1 ^2 



< rt>\p>)\ub-\B-{p) >= [(p + p% - ^ ™*\ \F*~*\f) + — ^" ^"'(g 2 ) (9) 

Other matrix elements are related to the above matrix elements using the following isopin 
relations. 

< n~(p)\du-\0 >=V2< 7T°\UU_\0 >= -V2< 7r°|cM-|0 > 

< K~(p)\su^\0 >=- < K°\sd-\0 >=- < K°\ds.\0 >, < ri {,) \uu-\0 >=< r] {,) \dd-\0 > 

< ir-(p')\db_\B-(p B ) >=< n + \ub_\B° >=V2< 7i°\db_\B° >= V2 < TT°\ub-\B- > 

< K-(j/)\sb-\B-(p B ) >=< K°\sb^\B° >, < r} {,) \ub_\B- >= - < r} {,) \db_\B° > (10) 

We choose the nemerical values of the relevant decay constants, f w , f K as follows: 

f n = 131 MeV, f K = 160 MeV. (11) 

The values of form factors at q 2 = ml is required for the calculation of decay rate. As 
the final states involve only light hadrons we can safely neglect the q 2 dependence of form 



5 



factors. Hence we assume that Fo^(m\) = _Fo,i(0). The values used for the form factors are 
obtained using the BSW model 



Fg{* K = 0.33, and Fg^ = 0.33 (12) 
In the case with nonleptonic B decays into the final state involving 77 or rj , we use two- 



angle mixing formalism developed by Leutwyler |T[] . This formalism is phenomenologically 
adequate to explain the various experimental data, the decay width r(?7 — > 27), T(rf — > 27) 
and the ratio r(J/ip — > rj'j)/T(J/i/j — > 777), as recently discussed in Refs. |21|J22]] . In this 
formalism, the 77 and 7/ states are defined by the mixture of the octet and singlet states with 
different mixing angles such as 

\r] >= cos^sl?^ > — sin^olr/o >, \r)' >= sm9 8 \r] 8 > + cos9 \r] > • (13) 

Then, the decay constants f s {l) are obtained through the two mixing angles and f 8 , fo- 



fu _ 
Jtj — 

f = 



-2^L cos 8 - -fa sin 9 , # = -2^ sin 9 8 + ^ cos 9 Q . 1 } 



And the relevant form factors for B — > 77 and B —* rf are: 



rpB^rr cos gs sin ftp rp^-* 1 !' pB- >7r sinflg j_ cos ftp /-i e\ 

^0,1 - ^0,1 ,/fi ,/a > -^0,1 - -^0,1 ,/fi ~+" ,/a • l iC V 



For numerical analysis, we choose the following values: 



9 8 = -22.2°, 6 = -9.1°, j- = 1.28, j- = 1.20 (16) 

Jit Jit 

which are given by fitting of data on the decay width T(rj — > 27), T(r]' — > 27) and the ratio 
T(J/ip — > 7/7) /T(J/ip — > 777) |p2| . The parameters j^, quantify the contributions from 
the decay 6 — > s(cc) — > 5(77,77'). Their magnitude and the sign were estimated using the 
QCD anomaly method in Ref. 0, and we use their values in this article : = —0.9 MeV 
and = -2.3 MeV for m c = 1.5 GeV. 

For numerical calculation of the transition matrix elements, we also need numerical values 



of CKM elements and quark masses. We use the Wolfenstein parameterization |23] of CKM 
matrix element. Two parameters A, A are well determined using \V c b\ through the fitting 
of the B — ► D*lvi decay spectrum and \V US \ through the K — ► nev and hyperon decays : 
A = 0.81 ± 0.06 and A = sin# c = 0.0025 ± 0.0018. In this paper, we choose the central 
value of these parameters. Other two parameters in CKM matrix elements are constrained 
by CKM unitarity fitting (HJ as 0.025 < 77 < 0.52 and -0.25 < p < 0.35 ( 95 % C.L. ). 
However the lower bound on the mass mixing ratio AM S / AM^ and the experimental value 
of #1 = B(B°(B°) -> n ± K^)/B(B ± -> vr ± ir ) = 0.65 ± 0.40 disfavor negative p region. 
Using the relation, y/ p 2 + rj 2 = \V u b\/X\Vcb\, we choose the three typical values of p and rj 
satisfying the above constraint following ref. [f|, 

(p = 0.05, 77 = 0.36), (p = 0.30, 77 = 0.42), (p = 0.00, 77 = 0.22). 

These correspond |Kd>|/|Vy = 0.08,0.11, and 0.05 representing the center value and upper 
and lower limit of the experimental value ( 90 % C.L. ) of the ratio. 
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In the calculation of Wilson coefficients, we use the internal quark masses as constituent 
quark masses: m b = 4.88 GeV, m c = 1.5 GeV, m s = 0.5 GeV, m d = m u = 0.2 GeV. 
The mass terms in the matrix elements of several decay modes come form the equation of 
motion. Hence they are current quark masses rather than the constituent quark masses and 
the values of them are given by m& = 4.88 GeV, m c = 1.3 GeV, m s = 122 MeV, = 7.6 
MeV, m u = 4.2 MeV at the renormalization scale \x = 2.5 GeV. 



III. BRANCHING RATIOS 

Armed with the strategies described in the previous section, it is straightforward to calcu- 
late the amplitudes and the branching ratios of B meson decays into two light pseudoscalar 
mesons. The full amplitudes are given in the Appendix, and we will consider the numerical 
results only in this section. We present the results for the averaged branching ratio, since 
CLEO has measured the averaged one : for example, 

B(B ± -> tjtt*) = - [B(B~ -> r]7T-) + B{B + -> r]n + )] . (17) 

In Table 1, we present the numerical results using the typical parameters : £ = 1/3, 
/(») = -0.9 MeV, ffi = -2.3 MeV ( m c = 1.5GeV ), p = 0.05, r) = 0.36. In Figs. 1, 2 and 
3, we plot the branching ratio as a function of nonfactorization parameter £ and using three 
different set of CKM angle p, rj. 

The additional effects of EW penguin diagrams in B — > tttt, ttK, KK modes, except 
Bd — > 7T 7r°, 7t°K° and B ± —>■ ir°K ± modes, are generally negligible and our predictions on 
the decay rates including EW penguin on the several modes are similar to that of Ref. [|J . 
and consistent on the recent CLEO data In fi ± — ► K^rj', the theoretical prediction 
has about 2a deviation from the CLEO data. In Bd — > 7r°7r°, 7r°K° modes, the EW penguin 
effects decrease the decay rate about 33% and 25% for £ = 1/3, respectively. In B ± — > ir°K ± 
mode, the decay rate is increased by the effects of EW penguin diagrams about 29% for 
£ = 1/3. The decay rate of B ± — ► K°K ± is same as that of Bd — * K°K° because of the 
isospin symmetry. 

In B — > Trr](\ Kt)(\ rj^'r]^ modes except B ± — > K ± r/,Bd — * K°r],r]7] modes, the EW 
penguin effects are also negligible and give at most C(few %) corrections in decay rates. For 
Bd — > T]T] modes, the effects increase the decay rate about 20%. In the B ± — > K ± r], Bd — > 
K°r], the decay rates are decreased by the effects about 39% and 36% respectively. In 
Bd — ► 7T 7r°, rf^rf' modes, the branching ratio plot has minimum value between £ = 0.2 
and £ = 0.3 and its values become very small in this region and very sensitive to the actual 
values of £, so that we cannot trust our predictions too much. 

We also compare the decay rates with and without the b —>■ s(cc) —>■ s(r), r}') contributions. 
These corrections does not change the decay rates significantly. Our results are quite different 



from those obtained in Ref. [25], where b — > s(cc) — * s(r},r)') effect was interpreted as the 
intrinsic charm contents in rj, rf mesons. The authors of Ref. [f25j estimated the numerical 
value of /i c V/n to be 0(1), and predicted that the branching ratios are about 0.194 x 
10~ 5 (5.83 x 10~ 5 ) for the decay modes B ± -> K^-qirf), and 0.027 x 10^(5.73 x 10" 5 ) for the 
decays modes B d — > K i](r]'), which are substantially larger than our predictions. However, 
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such large / ^ are not compatible with other theoretical/experimental considerations |5,22 
and cannot be taken too seriously. 



IV. CP ASYMMETRIES 



We also consider the direct CP violating rate asymmetry in the charged B meson decay 
and the time integrated CP asymmetry in neutral B meson decay. In the inclusive charmless 
B decays, the necessary strong phase differences are obtained from the absorptive part of 
the penguin diagram |19| , p6| , |27| . Such method is applied to exclusive nonleptonic decays 
of B meson decays into two light pseudoscalar mesons ||28|| . It had been usually assumed 



that the final state interactions in nonleptonic B decays may be negligibly small, because 
the mass of B meson is far above the usual resonance region. However such effects become 
important in some cases. For example, the soft rescattering effects |2D|-^TJ might affect the 
method to constrain and determine the CKM angle 7 using B — > irK modes p^ . 

In this article, we neglect the final state interactions and calculate the typical values and 
various parameter dependences of CP-asymmetries. 

The CP- violating rate asymmetry is defined as 



acp 



T(B- -> f) - T(B+ -> /) 

r(s--/) + r(fl+-/) 



(18) 



where / is the chrege conjugated state of the final state /. 



The time integrated CP-asymmetry of neutral B d meson |33 

1 



is 



a C p(B d -> /) 



1 + xj 



A% r P (B d - /) + A™p~ (B d - /) 



with 



AZ ix ~ ind {B d -> /) 



-2x d Imjf 

1 + IOI 2 ' 



and 



x d 



Am*, 



B d 



0.71, 







q A(B° -> f) 
P A(B° -> /)' 



(19) 



(20) 



(21) 



where / is CP eigenstates and Am denotes the mass difference of B° and B° mesons. The 
ratio q/p represents the CKM matrix elements ratio contributing to the B° — B° mixing. If 
we consider B meson decay into the final states with K°(K°) meson, the parameter £y has 
additional factor (q/p)K for the K° — K° mixing effects. 

We present the typical values of CP asymmetry of several modes and estimate the 
EW penguin effects and effects in Table 2. We also consider the effects of different 
CKM angles p, 77. In the Figures 4,5 and 6, the CP asymmetry is given by the function of 
nonfactorization parameter £ with three different CKM angles p, r]. In Ref pBf , the authors 
present the results of CP asymmetry as a function of q 2 . We will fix this value as ml/2 in 
this analysis. The EW penguin effects and the b — > d(s)[cc —> gg — » r/^] type corrections 
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on CP violation are generally small. The two-mixing scheme also gives generally small 
corrections to the values of the CP asymmetry compared with the one-mixing formalism for 
7] — rf system. 

In the CP asymmetry in charged B decays, The nonfactorization parameter £ dependence 
is generally mild. The CP asymmetry in B — > tt^tt , K 7t° is very small, at most 0.2% 
in magnitude. For B ± — > K^tc , K ± t] mode, the magnitude of CP asymmetry is less than 
about 10% in the most parameter sets. In B ± — > K ± K°, it^rf \ and K ±r q modes, the range 
of the magnitude of CP asymmetry is 10 ~ 30%. In the CP asymmetry in neutral B meson 
decays, B^ — > n°ir and B^ — > rf ^rf -* modes give very steep £ dependence in 0.2 ~ 0.4 region 
of £ region. The sign of CP asymmetry is changed in the region near £ ~ 0.3. Hence the 
typical values for such modes in table 2 ( which values are given in £ = 1/3 ) should be 
largely changed in small shift of £. In other mode, the CP asymmetry change slowly in 
varying f . 



V. POSSIBLE EFFECTS THROUGH ENHANCEMENTS IN B — > SG 

For the last several years, there have been some speculations about the enhanced b — > sg 



by several authors JT8[ in order to resolve some discrepancies in the data and theoretical 
expectations in inclusive B decays. The CLEO and ARGUS collaboration @,^] have 
measured the semileptonic branching ratio: 

B e sl = 10.23 ± 0.39% (22) 



And the CLEO 1.5, CLEO II and ARGUS data [31] give 



n e c xp = 1.15 ± 0.05 (23) 

for the average number of charm(anti-) quarks per B + / B°-decaj. The result ( f22"D is consid- 
erably smaller than the theoretical prediction in the parton model. It has been found that 
charm mass corrections to T(b — > ccs) are large and can reduce the theoretical prediction 
for B SL [H]: 

B SL = (H.7 ±1.4 ±1.0)% (24) 



Then, there is no spectacular discrepancy between (|22|) and (|24j). However, as we lower 
the theoretical prediction for B$l by increasing T(b —>■ ccs), we simultaneously increase the 
prediction for n c . The theoretical prediction for n c obtained from eq.(B3) reads 



)!,. 



1.34 =f 0.06 (25) 

The discrepancy between ( |2"3"| ) and ( ^5] ) constitutes the "missing charm puzzle". 

Furthermore CLEO collaboration reported rather large branching ratio of B — > r]'X s 
decay mode p7fl : 



B(B -> rj'Xs) = (6.2 ± 1.6 ± 1.3) x 10~ 4 for 2.0 < P n > < 2.7GeV (26) 
These anomalies may be solved through an enhancement of the chromomagnetic dipole 



coefficient C g by new physics. For example, Lenz et.al |38] showed that |C s (Mjy)| must be 
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enhanced by a factor of 9 to 16 in order to explain the observed charm deficit if the CKM 
structure of new physics contribution is the same as in the Standard Model. 

It is straightforward to include such effects of the enhanced chromomagnetic moment 
b — > sg on the nonleptonic two body decays of B mesons. This is accomplished by a simple 
replacement of 

Cg — > Cg = Cg^SM + C 'g <n eW ■ (27) 

Here, to see the possible effects on the nonleptonic two body decays of B meson, we take 
two possibilities : 

C g (2.5GeV) = +5 C g ,sM(2.5GeV), and - 5 C giS M(2.5GeV) (28) 

In Fig. 7 and Fig. 8, we show the predictions for B — > ttK and B — > r\K with 
Cg^new -^-^Cg t SM at pi — 2.5 GeV. The agreement between the data and the factorization 
predictions is generally degraded when one includes the enhanced b — > sg. For the positive 
,SMi all the predictions become worse because of destructive interference 
between the SM amplitude and the enhanced b ^ sg. For the negative C g ^ new = —5C 9t sM, 
one can improve the branching ratio for B — > 77 K , only by paying a price to the worse 
predictions for B — > 7r.fr. All these observations are based on the generalized factorization 
assumption. In such an approximation, there is a tendency that the enhanced b — > sg makes 
worse the agreements between the data and predictions. It has to be kept in mind that the 
enhanced b — > sg scenario should be also tested in the exclusive B decays in a better way 
than considered in this work, if possible. 



VI. CONCLUSIONS AND DISCUSSIONS 

In this work, we considered the branching ratios and CP asymmetries in B meson decays 
into two light pseudoscalar mesons, and the possible effects of the enhanced b — > sg vertex 
suggested as a solution to the semileptonic branching ratio problem and the missing charm 
puzzle in B meson decays. The typical branching ratios in various parameter set are given in 
Table 1. Their £ parameter dependences are given in Figs. 1, 2 and 3. In B — > tttt, ttK, KK 
modes, we included the EW penguin effects as well as gluonic penguin effects. Especially 
in B d — > -K -K ,7c°K , K°r] modes and B ± — > tc ^, K ± r\ modes, EW penguin effects give 
important contributions to the branching ratio: about +33%, +25%, —36%, +29%, —38% 
variation respectively. In B& — > 7r°7T ,?/ ) modes, the branching ratio plot has minimum 
value between £ = 0.2 and £ = 0.3 and its values become very small in this region and very 
different from the values outside such the region. 

The EW penguin effects and the b — > d(s)[cc — > gg — > ?/')] type corrections on CP 
violation generally give small contributions to CP asymmetry. The typical values of the CP 
asymmetry of B meson decays are given in Table 2. In Figs. 4, 5 and 6, the £ dependences 
of CP asymmetry in several modes are given in three different set of CKM angle p, r]. For 
a fixed £, the important corrections of CP asymmetry are given by varying the CKM angle 
itself. The £ dependece of CP asymmetry in charged B decays are rather mild. In the 
neutral B meson decay modes, CP asymmetry of B d — > 7r°7r° and B d — > rj^rj^ modes 
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change very sharply in the region 0.2 < £ < 0.4. In other modes, the CP asymmetry 
changes slowly when one varies £. 

Lastly, we made an observation that the enhanced b — > sg tends to degrade the agreement 
between the data and the factorization predictions for B meson decays into two light mesons, 
which seems to disfavor the enhanced b — > sg scenario. The enhanced b — > sg idea was put 
forward within the inclusive B decays ||18|| , which can be studied with much less theoretical 
uncertainties compared to the exclusive cases. But it is certainly true that the enhanced 
b —>■ sg will affect the individual exclusive B decay modes as well. It would be very welcome 
to study the effect of this enhanced b — > sg on the individual exclusive B decays in a way 
more reliable than the generalized factorization method employed in the present work. 



VII. APPENDIX 



In this appendix, we present the complete matrix elements including the EW penguin 
effects and effects. 



B — ► 7T 7T 
G 

72 I 



M = S lv uh V: d ( ai + a 2 ) - V tb V; d 



3 , 2mla 8 

- a 9 - a 7 + a w ) + ■ 

2 (m u + m d )(m b - m u ) 



mla 8 



2m d (m b - m u ) 



< 7r \du_\0 >< ir"\ub-\B > 



(29) 



where 



< 7r-|du_|0 >< K°\ub-\B- >= % ^={m% - ml)F^(ml] 

v2 



(30) 



B° 



IX IX 



O4 + flio + 



(m u + m d ){m b - m l 



x < ix \du-\Q >< ix + \ub_\B° > 



(31) 



where 



< ir~\du-\0 >< ix + \ub-\B° >= i f n (m 2 B - ml)F^ n (m 2 n ) 



(32) 



IX IX 



M = ^{v ub V: d a 2 + V tb V t * d 



a 4 + -{a 7 - a 9 , 



1 q m^(2a 6 - a 8 ) 
2 2m d (m b - m d ) 



x < ix°\uu\0 >< ix°\db\B° > 



(33) 



where 



< tx°\uu_\0 >< ix°\db_\B° >= 1 —{m\ - ml)F^ n (ml). 



(34) 
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B~ -> K°7T 

V2 



M = —^V tb V t * s 



a - -a + raj-( 2a 6 - g 8 ) 

2 (m d + m s )(m b - m d ) 



< K u \sd.\0 >< n-\db-\B- > 

(35) 



where 



< K u \sd_\0 >< ir-\db_\B- >= -i f K (m 2 B -ml)F^ n (m 2 K ) 



(36) 



B~ -> ^-7T° 



O4 + Oio + 



2m 2 K (a 6 + a 8 ) 



+ -(a 9 - a 7 )R K - n o 



(m u + m s )(m b - m u ) 
I < fT|sw_|0 >< 7T°|u6_|B- > (37) 



where 



<7T°| 


uu_\0 >< K~\sb- 


B~ > 


< K~ 


\su-\0 >< 7r°\ub- 


B- > 



Rk-tt° = 



< K-\su_\0 >< 7v°\ub_\B- >=i ^=(m 2 B -ml)F^ n (m 2 



V2' 



K J 



(38) 



M=^iv uh V: s a 1 ~V tb V t * s 



a 4 + a 10 + 



V2 I 

x < tf-|su_|0 >< n + \ub_\B° > 



2m 2 K (a 6 + a 8 ) 
[m u + m s )(m b - m u ) 



(39) 



where 



< AT-|su_|0 >< it + \ub_\B° >= i f K (m 2 B -ml)F^{m 2 K ). 



(40) 



B° -> i?°7T 



a - -a + m ^( 2 ° 6 ~ a s) 

2 (m d + m s )(m b - m d ) 



< K°\sd_\0 >< ir°\db_\B° > 



(41) 



where 



R K o n o 



< tt°|mm_|0 >< ir°|s6_| J B° > _Um%- m\ F^ K {m\) 

< K°\sd„\0 >< 7T°\db.\B > ~ ~ fa m 2 B -m 2 n F B ^(m 2 K ) ' 



< X°|sd_|0 >< 7i°\db_\B° >= -i ^=(m 2 B - m 2 K )F*^(m 2 K ). 

a/2 



(42) 
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B~ 



K°K~ 



M = 



—V tb V* ja 4 - -aio + m K( 2a e ~ a s) 
V2 1 U \ 4 2 10 (m d + m s )(m b -m s ) 

x < K°\ds_\0 >< K~\sb-\B- > 



where 



< K°\ds_\0 >< K~\sb-\B- >= -i f K (m 2 B - m\)F^ K [m\) 



M = 



G, 



m 2 K (2a 6 - a 8 ) 



1 

2 (m d + m s ) (ra 6 - m s ) 
x < fT°|Js_|0 >< #°|s&_|fl > 



where 



< K°\ds_\0 >< K°\sb^\B° >= -i f K (m 2 B - m 2 K )F^ K (m 2 K ). 
B~ — > ix~r\ 
M = 

-v tb v; d 



Gf { 

yubV* d ( a 2 + ai R n - V ) + V cb V* d a 2 -j^ 



2(a 3 - a 5 ) + a 4 + 



m 2 (2a 6 -a 8 ) ( ft\ 1 



+ < a 4 + aio + 



2m^(a 6 + a 8 ) 
(m u + m d )(m b - m„) 



R 



TT~ri 



+ |a 3 - a 5 + -(a 7 - a 9 )| + {a 3 - a 5 - a 7 + a 9 } | 



x < tt |d6_| B ><ri\uu„\0> 



where 



R-n 



_ < 7T~ 


uu-\0 >< r]\db_ 


B~> _U m l~ m l 


57 < 77|tm_|0 >< 7T-\db- 


B~ > f% m 2 B - m 2 n 



< 7r-\db_\ B- >< V \uu_\0 >= i /« (m| - m\) F^{m 2 K ) 



B° -> TT /? 
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M 



^§ ^V ub V* d a 2 (1 + itVg + V cb V* d a 2 ^ 



-v tb v; d 



of \> m^(2a 6 -g 8 ) ( /«\ 1 



- { a 4 + (a 6 - -a 8 ) 



m: 



3 1 I 



r i l / s / (c) l) 

+ ja 3 - a 5 + -(a 7 - a 9 )| + {a 3 - a 5 - a 7 + a 9 )} > 
x < ri\uu_\0 >< 7T°\db-\B° > 



(49) 



where 



< tt°|mm_|0 >< r]\db_\B° > 



U ™b 



n v < r]\uu_\0 >< ir°\db_\B° > ~ /« m\-m\ F^{m^ 
< ri\uu-\0 >< n°\db-\B >= i (m\ - m\) F^{m\). 



(50) 



The matrix elements of the B — > 7r 77', 5° — » 7r°r/ modes might be obtained by 
replacing 77 with rj' in — > n^i], B° — > 7r°r? modes. 



M = % { VubKs ( ° 2 + + VcbV * s ^ " Ws 



2(a 3 - a 5 ) - 



m^(2a 6 - a 8 ) 
2m s (m b - m s ) 



z, ^ v ;i s < IIL u)\' IL b IIL u) ) 

( m\{2a % -a 8 ) 1 1 /» 



+ < a 3 + a 4 — a 5 + 



2 m ,( mi - m ,) + 2 (ai - a »- ai »»|7» 



+ (a 3 - a 5 - a 7 + a 9 ) 



f(c) 
j »7 



< r)\uu-\0 >< K \sb-\B > 



(51) 



where 



__ < K~ 


su-\0 >< r)\ub- 


B- > _ f K m%- m\ 


v < rj\uu_\0 >< K~\sb_ 


B~ > /* m%- m\ 



< r]\uu^\0 >< K-\db-\B~ >= i f% (m 2 B - m\) F^ K {m\). 
B° -> K°r] 



M 



^\v ub V: s a 2 + V cb V: s a 2 ^ 



- v tb v* 



2(a 3 - a 5 ) - 



m^(2a 6 - a 8 ) 
2m s (m b - m s ) 



(52) 



-(a 7 -a 9 ) 
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+ < a 4 + 



m 2 K (2a 6 - 



+ < a 3 + a 4 — a 5 



(m s + m d )(m b -m d ) 2 

m^(2a 6 -a 8 ) 1 



f s 



2m s (m b -m s ) ' 2 ( ° 7 " ° 9 " aio) | # 



+ (a 3 — a 5 — a 7 + a 9 ) 



f(c) 



< r]|MM„|0 >< A"°|s&_|.B > 



(53) 



where 



= < if°|sd_|0 >< r]\db_\B° > _ f K m\-m\ F^\m 2 K ) 



< r]\uu_\0 >< K°\sb_\B° > /« m| - F S " 
< 77|uu_|0 >< >= i /; (m%-m 2 K ) F M (mJ). 



(54) 



The matrix elements of the B — > X 77', S° — > fT ^' modes might be obtained by 
replacing 77 with 77' in I? - — > -B — > -fT°77 modes. 



5° -> 7777 

M = 



^\v ub V: d a 2 + V cb V; d a 2 f -^ 



(55) 



id 



2(03 " + ° 4 + 2777,(777,-77,,) ~ f) " 2^ " ^ + ^ 



1 



f 



f \ a 3 - a 5 + -(a 7 - a 9 ) j j| + (a 3 - a 5 - a 7 + a 9 ) 
x < 77|mm_|0 >< T^l^li? > 



, ■ / »7 



c) 



Jr; J 



(56) 



where 



< 7?|7m_|0 >< r]\db^\B° >= -i f% (m 2 B - m\) F^ v {m 2 v ) 



(57) 



The matrix elements of the B° — > rj'i]' might be obtained by replacing 77 with 77' in 
B° — > 7777 modes. 



B° — > 7777' 
M 



Gp { f (c) 

m 2 {2a 6 



(58) 



-v tb v t 



id 



2(a 3 - a 5 ) + a 4 + 



2m s (m b - m s ) \ f* J 2 
I 1 f s f( c ) 



(a 7 - a 9 + a w ) 



+ |a 3 - as + 2 (°7 - a 9) J + (a 3 - a 5 - a 7 + a 9 ) 
+ (77 -> 77') 



< r]\uu^\0 >< ri'\db^\B° > 

(59) 
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where 



< ti\uu_\0 >< rf\db-\B° > = -i f% (m| - rrfy) F W (mJ), 

< rf \uu-\0 ><7i\db_\B°> = -i ffi (m| - ml) F^ v {ml). (60) 
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TABLES 



TABLE I. Combined branching ratios in unit of 1(T 5 . 'QCD' and 'EW present the QCD 
penguin and EW penguin effects respectively. 'Full' includes the b — > s[cc — > gg — ► rjrf] effects : 
£ = 1/3, 4 c) = -0.9MeL, ffi = -2.3MeV ( m c = 1.5GeV ),p = 0.05, r? = 0.36 



Decay Mode 


Tree 


Iree+QUD 


Tree+QCD+EW 


i ' , , . 

_hxp. 


_D — > 7T 7T 


0.51 




U.Ol 


0.51 








£> — > 7T A 







l.oo 


1.61 






O Q + l-1 i no i no 

l.6_i q ± 0.6 ± U.Z 


d — > 7T n 


0.038 




U.o4 


1.08 






<1.D 


£> — > A A 







U.Uoz 


0.081 




<2.1 


£>(i — > 7T 7T^ 


0.94 




u.yo 


0.93 






<1.0 


D „0„0 


0.11 x 10" 


-2 


n noi 

U.Uzl 


0.014 




n no 


±5^ — > 7T A ' 


0.071 




l.oo 


1.71 






-i r+0.5 in i 1 n 1 

1.0_q 4 ± U.l zt U.l 


B d -> 7T U A U 


0.53 x 10" 


-4 


0.83 


0.62 






<4.1 


_, ^0 
£>(i — > A A 







U.UoZ 


0.081 






Decay Mode 


Tree 


Tree+QCD 


Tree+QCD+EW 


Full 


Exp. 


B ± -> 71^7? 


0.23 




0.25 


0.26 




0.26 


<1.5 




0.16 




0.18 


0.17 




0.18 


<3.1 


£d A 


0.18 x 10" 


-4 


0.026 


0.026 




0.026 


<0.8 


B d - A' 


0.71 x 10- 


-5 


0.011 


0.79 x 10~ 2 




0.81 x 10~ 2 


<1.1 


£± _^ 


0.017 




0.18 


0.11 




0.11 


<1.4 




0.012 




2.15 


2.05 




2.08 


6.5±|;f ±0.9 




0.37 x 10" 


-4 


0.15 


0.096 




0.094 


<3.3 


5 d - #V 


0.24 x 10' 


-4 


2.16 


2.04 




2.07 


4.7± 2 ;£ ± 0.9 


Bd -> W 


0.32 x 10" 


-3 


0.83 x 10" 2 


0.010 




0.010 


<1.8 


-> 777/ 


0.44 x 10" 


-3 


0.78 x 10~ 2 


0.83 x 10~ 2 




0.86 x 10~ 2 


<2.7 


£ d - 7/7/ 


0.15 x 10" 


-3 


0.17 x 10~ 3 


0.14 x 10~ 3 




0.15 x 10~ 3 


<4.7 
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TABLE II. CP asymmetries in %. 'QCD' and 'EW present the QCD penguin and 
EW penguin effects respectively. 'Full' includes the b — > s[cc — ► gg — ► 7777'] effects. 
'Full(one-mixing)' presents the one- mixing scheme for 77 — 77' mixing : CKM phases are 
p = 0.05, 7] = 0.36(p = 0.30, 77 = 0.42; p = 0, 77 = 0.22), £ = 1/3 



Decay Mode 


QCD 


QCD+EW 


D± ± 

B^ — ► 7T 7T 
±5 — > 7T K 


( 0, 0) 
1.6 ( 1.8, 1.0) 




-0.05(-0.03,-0.09) 
1.6 ( 1.8, 1.0) 




r>+ jvH- 
— > 7T U if ± 


8.4 ( 12.9, 5.0) 




6.7 ( 10.0, 4.0) 




±? — > A U A =C 


— 12. 2(— 20.5, -7.3) 




-12.3(-20.6,-7.3) 




r> -(- zr 
— > 7r i 7T+ 


5.1 ( 21.4, 22.8) 




5.4 ( 21.5, 23.3) 




D 

B d -> 7r u vr u 


— 10. 5(— 18.8, —6.5) 




-13.9(-24.4,-8.6) 




d Or/ 


31.1 ( 41.1 19.8) 




31.1 ( 41.1, 19.8) 






15.2 ( 26.0, 9.1) 




15.2 ( 26.1, 9.1) 




Decay Mode 


QCD 


QCD+EW 


Full 


Full(one- mixing) 


S± 71^77 


-18.2(-9.6,-23.5) 


-18.0(-9.6,-22.4) 


-18.0(-9.6,-22.4) 


-18.6(-9.9, -23.1) 


B ± TT*??' 


-17.9(-9.4,-25.9) 


-18.0(-9.4,-26.3) 


-17.9(-9.4,-26.1) 


-17.7(-9.3, -25.9) 


— v ft^n 


_fi 7(—^ "i —A R\ 

U. 1 l tJ.O, ^i.O 1 


-10.9(-7.9,-8.2) 


-ll.l(-8.0,-8.4) 


-13.4(-7.5,-14.5) 




4.5 ( 5.8, 2.7) 


4.7 ( 6.1, 2.8) 


4.6 ( 6.0, 2.8) 


4.8 ( 6.3, 2.9) 


B d - A 


18.8(31.4,11.3) 


18.8(31.4,11.3) 


18.9(31.6,11.3) 


19.2(32.0,11.5) 


^ - A' 


23.5(38.4, 14.2) 


27.9(44.5,16.8) 


28.1(44.8,17.0) 


28.1(44.7,17.0) 


B d - ^577 


31.8 ( 41.7, 20.3) 


32.0 ( 41.9, 20.5) 


32.1 ( 41.9, 20.5) 


34.3 ( 43.6, 22.1) 


£ d - 


29.9 ( 40.1, 19.0) 


29.9 ( 40.0, 18.9) 


29.9 ( 40.1, 18.9) 


29.9 ( 40.0, 18.9) 


Brf -> 7777 


37.3 ( 52.4, 22.9) 


33.1 ( 48.1, 20.2) 


33.1 ( 48.1, 20.2) 


33.5 ( 48.6, 20.4) 


Bd -»■ W 


45.5 ( 59.6, 28.5) 


43.9 ( 58.3, 27.4) 


43.8 ( 58.1, 27.4) 


43.9 ( 58.3, 27.4) 




57.7 ( 67.8, 37.8) 


63.1 ( 70.6, 42.5) 


62.4 ( 69.9, 41.9) 


62.5 ( 69.8, 42.0) 
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FIGURES 
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1 
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upper boundary of Exp. 



lower boundary of Exp. 
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(d) B" -> K° 



3.1 I upper limit ( 90 % C.L. ) 



r T T--i--i--i---i---i---i---i---i"-r"r" i i i i i i i 

0.2 0.4 0.6 0.8 1 



(0 



7T 7] 



p 



6 
A 
2 




upper boundary of Exp. 



lower boundary of Exp. 



2~L 
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FIG. 1. Combined branching ratio of the charged B meson decays as a function of nonfac- 
torization parameter £: p = 0.05, r/ = 0.36 for real line, p = 0.30, r/ = 0.42 for dashed line and 
p = 0, 7] = 0.22 for dotted line. 
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FIG. 2. Combined branching ratio of the neutral B meson decays as a function of nonfac- 
torization parameter £: p = 0.05, 77 = 0.36 for real line, p = 0.30, r] = 0.42 for dashed line and 
p = 0, 77 = 0.22 for dotted line. 
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FIG. 3. Combined branching ratio of the neutral B meson decays as a function of nonfac- 
torization parameter £: p = 0.05, 77 = 0.36 for real line, p = 0.30, r] = 0.42 for dashed line and 
p = 0, 77 = 0.22 for dotted line. 
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FIG. 4. The CP asymmetry of the charged B meson decays as a function of nonfactorization 
parameter £: p = 0.05, r] = 0.36 for real line, p = 0.30, r] = 0.42 for dashed line and p = 0, r] = 0.22 
for dotted line. 
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FIG. 5. The CP asymmetry of the neutral B meson decays as a function of nonfactorization 
parameter £: p = 0.05, r] = 0.36 for real line, p = 0.30, r] = 0.42 for dashed line and p = 0, r] = 0.22 
for dotted line. 
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FIG. 6. The CP asymmetry of the neutral B meson decays as a function of nonfactorization 
parameter £: p = 0.05, r] = 0.36 for real line, p = 0.30, r] = 0.42 for dashed line and p = 0, r] = 0.22 
for dotted line. 
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FIG. 7. Combined branching ratios of the B meson decays as functions of nonfactorization 
parameter £ in the presence of the enhanced b — > sg with C g ^ new = +5C gj sM at /x = 2.5 GeV : 
p = 0.05, 7] = 0.36 for real line, p = 0.30, r) = 0.42 for dashed line and p = 0, rj = 0.22 for dotted 
line. 
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FIG. 8. Combined branching ratios of the B meson decays as functions of nonfactorization 
parameter £ in the presence of the enhanced b — ► sg with C g ^ new = —5C gj sM at // = 2.5 GeV : 
p = 0.05, 7] = 0.36 for real line, p = 0.30, r) = 0.42 for dashed line and p = 0, rj = 0.22 for dotted 
line. 
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